Negative four-terminal resistance as a probe of crossed Andreev reflection 
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We report on the experimental investigation of electronic transport in superconductor-ferromagnet 
spin- valve structures. Our samples consist of two ferromagnetic iron leads forming planar tunnel 
contacts to a superconducting aluminum wire. At energies below the superconducting gap, we 
observe a negative four-probe resistance that can be explained by crossed Andreev reflection. 
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Quantum mechanics predicts a non-local correlation 
of spatially separated particles that is stronger than that 
allowed by classical theory. This correlation, known as 
entanglement, has been tested experimentally on pairs of 
photonsi and massive particles^, see alsoi An intrigu- 
ing question is whether a stream of spatially separated, 
entangled electron pairs can be created and manipulated 
in a solid-state environment, with possible future appli- 
cations in quantum information processing. A promising 
implementation of solid-state entanglers are heterostruc- 
tures involving a superconductor;^ as in a superconduc- 
tor the electrons come naturally in the form of Cooper 
pairs, which are in an entangled singlet state. The chal- 
lenge is to create structures where the two electrons of a 
Cooper pair can be spatially separated without destroy- 
ing entanglement. If two normal- metal contacts are at- 
tached to a superconductor at a distance smaller than the 
coherence length of the superconductor, it has been pre- 
dicted that an electron injected into the superconductor 
may be transmitted into the second contact as a hole, 
thereby creating a Cooper pair in the superconductor, 
and leaving behind two entangled holes in the two nor- 
mal metal contacts^. There is a competition between 
this process, called crossed Andreev reflection (CAR), 
and elastic cotunneling (EC), where the incident electron 
gets transmitted to the second contact via a virtual state 
in the superconductor. Quantitative understanding of 
the contributions of CAR and competing processes like 
EC is a prerequisite to building a solid-state entangler 
based on superconductor hybrid structures. 

Crossed Andreev reflection has been studied ex- 
perimentally in superconductor /ferromagnet* and 
superconductor /normal- metal* structures. In a recent 
experiment^ on transport properties of superconduc- 
tor/ferromagnet spin- valve structures, we have observed 
spin-dependent transport in the superconductor at 
energies much smaller than the superconducting gap. 
The observed signal decayed on the length scale of the 
coherence length of the superconductor, and could be 
explained by theoretical predictions for the superposition 
of crossed Andreev reflection and elastic cotunneling. 
Here, we present an experimental study that allows us 
to discriminate the two processes, and also rule out 
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FIG. 1: (Color online) Schematic view of the experiment. 
Two ferromagnetic leads form tunnel contacts to a supercon- 
ducting bar. Contact A is used to inject a current, and the 
voltage at contact B is measured referenced to the drain volt- 
age probe. The magnetization alignment of the ferromagnets 
is either parallel or antiparallel. 



sequential tunneling or non-equilibrium effects as a 
possible cause of the observed signals. 

The schematic view of the experiment is shown in Fig. 
^ Two ferromagnetic wires A and B are attached via 
tunnel contacts to a superconducting bar. Contact A 
(source) is used to inject a current into the superconduc- 
tor, which then flows along the bar. The (Ohmic) drain 
contact, which lies inside the current path, is used as a 
reference of the chemical potential of the superconductor. 
Contact B is used to measure voltage relative to the drain 
contact. There are several possible mechanisms for the 
observation of a finite voltage at contact B despite the 
fact that voltage is measured across a superconductor: 
EC (CAR) emits an electron (hole) from the supercon- 
ductor into contact B for an electron injected into contact 
A, thereby leading to a voltage that has equal (opposite) 
sign compared to the injector voltage. Both processes 
are elastic and coherent, and involve virtual quasiparti- 
cle states in the superconductor, i.e. they may occur at 
energies below the superconducting energy gap, and over 
the length scale of the coherence length £s of the super- 
conductor. Note that for CAR, the chemical potential of 
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FIG. 2: SEM image of the two contacts of sample II. 



the detector contact B may lie outside the voltage win- 
dow spanned by source and drain contact, which has been 
predicted previously for a slightly different setup^S 

In addition to coherent subgap processes, an electron 
may tunnel into an allowed state in the superconduc- 
tor and sequentially tunnel out at contact B. Sequen- 
tial tunneling, like elastic cotunneling, yields a voltage 
Vb of the same sign as Va- Sequential tunneling re- 
quires an electron to have enough energy to overcome the 
spectral gap, which can be supplied by the applied volt- 
age, thermal excitation or external noise. Alternatively, 
the gap may be smeared out by quasi-particle life-time 
broadening, the inverse proximity effect in the contact 
regions, or the magnetic stray fields of the electrodes. 
Sequential tunneling may be either elastic and coherent, 
or incoherent, generating a non-equilibrium charg o 1 ^ 12 
and spirit accumulation in the superconductor. It is 
known that non-equilibrium quasiparticlc populations in 
superconductors relax only slowly due to electron-phonon 
scattering^ and accordingly the decay length of the non- 
equilibrium can be quite large. From the data of our 
previous experiment^ we obtain the charge imbalance 
relaxation length Aq* « 3 /im and the normal-state spin 
diffusion length A s f w I /im. A 10 /im distance between 
source and drain contact was chosen, to ensure that any 
non-equilibrium quasi-particle populations which may be 
injected at contact A have relaxed at the drain reference 
contact. Therefore, the chemical potentials of both quasi- 
particle spin species (^qp ) and the Cooper pairs fi cp are 
the same, and the drain contact can be considered as a 
faithful measure of the equilibrium chemical potential of 
the superconductor. 

The experiments were performed on mesoscopic spin- 
valve structures made of ferromagnetic iron and super- 
conducting aluminum by means of e-beam lithography 
and shadow evaporation. The contact region of sam- 
ple II is shown in the SEM image in Fig. [21 In a first 
evaporation step, 15 nm of iron is evaporated from a di- 
rection almost normal to the plane of an oxidized silicon 
wafer. The purpose of this first iron layer is only to en- 
sure Ohmic contacts everywhere in the structure except 
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FIG. 3: (Color online) Temperature dependence of the de- 
tector voltage Vb at fixed injector current 7a = 11 nA for 
parallel (diamonds) and antiparallel (circles) magnetization 
alignment for sample II. Inset: Vb as a function of 7a in the 
antiparallel alignment at T — 25 mK. A marks the current 
that corresponds to Va = 200 fiV, i.e. the gap energy of 
aluminum. 



for the two tunnel contacts A and B. Then, 25 nm of alu- 
minum are evaporated under an oblique angle, forming 
the aluminum wire of 200 nm width seen in the upper 
part of the SEM image. The aluminum is then oxidized 
in situ admitting a pressure of 60 Pa of pure oxygen for 
about 1 min. Subsequently, a second layer of 20 nm of 
iron is evaporated from the opposite direction to form the 
two tunnel contacts A and B. After lift-off, the samples 
were bonded and mounted into a shielded box thermally 
anchored to the mixing chamber of a dilution refrigera- 
tor. All wire connections from the top of the cryostat 
to the sample chamber were filtered through discrete RC 
low-pass filters, lossy coaxial lines and a copper pow- 
der filter to eliminate noise and high-frequency thermal 
photons from the sample chamber. The experiments pre- 
sented here were carried out using a battery and a resistor 
network as DC current source, and a nanovoltmeter for 
voltage measurements. The low-bias results were con- 
firmed using an AC resistance bridge. Here, we present 
data from three samples I, II and III. 

Figure |21 shows the detector voltage Vb as a function 
of temperature for sample II. Here, the magnetization of 
injector and detector contact have been aligned parallel 
or antiparallel at T > T c by monitoring the normal-state 
spin-valve signal, and then subsequently the sample is 
cooled down at zero magnetic field. Above T c , the ob- 
served voltage is essentially due to the Ohmic resistance 
of the 10 /im long aluminum strip between contact B 
and the drain reference contact, plus the small spin ac- 
cumulation signal. At T c w 1.35 K, the Ohmic resistance 
disappears, and a small positive signal remains around 
1 K. At lower temperature, the detector voltage becomes 
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negative for both magnetization alignments, but with a 
much larger magnitude for antiparallel alignment. 

The dependence of the detector voltage Vb on bias cur- 
rent /a is shown in the inset of Fig. [3] At low bias current 
/a | < 30 nA, the slope is negative, turning into a posi- 
tive slope at higher bias. By comparison with the local 
current-voltage relation (IV) of the injector contact we 
see that the slope reversal corresponds to Va ~ 200 /iV, 
i.e. the superconducting energy gap of aluminum. Thus, 
we conclude that the positive slope corresponds to the on- 
set of the transmission of quasiparticles through allowed 
states above the gap. 

From the measured non-local voltage as a function of 
current, the non-local conductance can be calculated by 
using the conductance matrix 



Ia 
Ib 



Ga Gx 
Gx Gb 



V a 
Vb 



(1) 



which contains the local conductances Ga,b of the in- 
jector (A) and detector (B) contact, and the non-local 
conductance Gx- For voltage detection at B, we have 
Ib = 0, and the current due to the non-local conduc- 
tance, Ix — GxVa, is canceled by the backflow due to 
the local conductance, GbVb- By balancing the total 
currents, we assume that all charge carriers emitted into 
contact B by non-local processes relax to the equilibrium 
chemical potential of wire B before they tunnel back. In 
that Vb never exceeds a few microvolts, we can 

replace Gb by its low-bias value, and the energy depen- 
dence of Vb is given only by the dependence of Gx on 
Va- With Gx < G A , we find 



G x 



dI B 
dV A 



-Gi 



dV B 
'dV A 



(2) 



The resulting non-local differential conductance 
dlB/dVA as a function of injector voltage Va, together 
with the local differential conductances of injector and 
detector contact measured directly, is shown in Fig. 0] 

The local conductance of the detector exhibits a sup- 
pression at subgap energies, and peaks at the energy of 
the gap, A « 200/xeV, which at least qualitatively cor- 
responds to the expected BCS tunneling characteristics. 
For the injector, neither the subgap suppression nor the 
peaks are well resolved. For both injector and detector, 
the subgap conductance is much too large to be compat- 
ible with a simple BCS behavior, and has a downward 
curvature at low bias, which is not expected for cither 
BCS quasiparticle tunneling, or Andreev reflection. A 
strong enhancement of Andreev reflection at low bias in 
mesoscopic NS tunnel junctions is known to occur due 
to the constructive interference of time-reversed electron 
and hole trajectories on the normal metal side of the tun- 
nel junction (reflectionless tunneling). 15 While reflection- 
less tunneling exhibits the decrease of conductance with 
increasing bias in the subgap regime that we observe, it 
should not occur in our structures, as time-reversal sym- 
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FIG. 4: (Color online) Local differential conductances 
dlA/dV\ and cUb/cLVb (left scale), and non-local differen- 
tial conductance cUb/cIVa (right scale), of sample II at T ~ 
25 mK. cLIa/cLYa is doubled for clarity. 



metry is broken in the ferromagnetic electrodes. How- 
ever, it has been predicted^ that an interference en- 
hancement of Andreev reflection also occurs on the super- 
conducting side of the interface in restricted geometries. 
The formalism used to describe this enhancement is the 
same as used later to describe crossed Andreev reflection 
in diffusive metalsiiii^ 

The non-local conductance is positive (corresponding 
to crossed Andreev reflection) in the subgap regime, and 
becomes negative (corresponding to electron transmis- 
sion) as the injector voltage approaches the supercon- 
ducting gap. At the gap, a negative peak is observed, 
mimicking the BCS density of states peak seen in the 
local conductance. The overall magnitude of the non- 
local conductance is about a factor of 100 smaller than 
the local conductance. Using the model of twodimen- 
sional diffusioni£*i&, we find that a relative reduction of 
the non-local conductance compared to the local An- 
dreev conductance of two orders of magnitude is real- 
istic. However, using the full quantitative expression^, 
the local Andreev conductance should be of the order of 
G^Rn 1 /xS, where Gn is the normal-state tunnel con- 
ductance, and Ra the normal-state sheet resistance of 
the superconductor. This estimate is much less than the 
conductances actually observed. We suspect that pin- 
holes in the oxide barrier lead to an enhanced Andreev 
conductance. 

Figure shows the non-local conductance as a func- 
tion of injector voltage for three different samples in the 
antiparallel magnetization state. For samples I and II, 
the data are qualitatively the same, with a dominating 
positive signal at low bias, and a negative peak near the 
energy gap at about 200 /xV. However, sample III shows 
a different behavior. The non-local conductance is neg- 
ative at small bias, has a positive peak at about 80 /iV, 
and then decays without any clear feature near the gap. 
For sample III, the detector voltage Vb has the same sign 
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TABLE I: Characteristic parameters of our samples. Average 
transmission probability t of the tunnel contacts, elastic mean 
free path i e i, coherence length £, normalized contact distance 
d/£, and inverse diffusion time H/tb- 
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8 x 1(T 5 


13 
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FIG. 5: (Color online) Non-local differential conductance 
dls/dVA as a function of injector voltage Va for three dif- 
ferent samples. The data are scaled to a similar amplitude. 

as the injector voltage over the whole IV trace. 

The theoretical predictions for the non-local conduc- 
tance due to CAR and EC in the tunneling limifpi 9 . show 
that both contributions are of the same magnitude and 
opposite sign. Therefore, without spin selection, the non- 
local subgap conductance is expected to be zero. In- 
cluding spin selection a negative (positive) conductance 
of equal magnitude is expected for parallel (antiparal- 
lel) alignment. For increasing contact transparency, EC 
is predicted to give a larger absolute contribution than 
CAR. 20 Here, we observe a dominantly positive non-local 
conductance for samples I and II, and a negative signal 
for sample III. In Tabled we summarize the experimen- 
tal parameters of the three samples. As can be seen, 
the average contact transparency estimated from the 
normal-state tunnel resistance systematically decreases 
from sample I to III. A transition from dominating CAR 
to dominating EC with decreasing transparency is in- 
consistent with the predictions of Ref. |2(J. The elastic 
mean free path l c \ of the aluminum also systematically 
decreases from sample I to III, and thereby the coher- 
ence length £ decreases and the normalized contact dis- 



tance d/£ increases. The dependence of CAR and EC on 
contact distance is expected to be the same ) 17 ' 19 which 
means that except for an overall signal decrease, no qual- 
itative change is expected with increasing c£/£. Recently, 
it has been predicted^! that coupling to the electromag- 
netic environment (dynamical Coulomb blockade) may 
discriminate CAR and EC, and also explain a transition 
between dominating CAR and EC at a finite bias volt- 
age related to typical energies of environmental modes, 
either symmetric or antisymmetric with respect to the 
two contacts. As the elastic mean free path of the alu- 
minum may affect the environmental modes, we specu- 
late that the qualitatively different behavior of sample III 
compared to I and II may be related to the decrease in 
Z e i, probably in conjunction with the decrease in contact 
transparency. 

The non-local conductance of sample III is qualita- 
tively consistent with the observations made by Russo 
et al^ in a coplanar geometry, where a superconduct- 
ing layer was sandwiched between two non-magnetic con- 
tacts. Russo et al. report a positive differential voltage 
(corresponding to a negative differential conductance) at 
low bias, which reverses sign at an energy scale well be- 
low the superconducting gap, and decays towards higher 
energies. The energy scale of the sign reversal was found 
to be the Thouless energy £rh = ft/ T D, where td is the 
diffusion time corresponding to the thickness of the su- 
perconducting film (which coincides with the distance be- 
tween the two normal metal contacts in a coplanar ge- 
ometry). For comparison, we have calculated the energy 
scale Ti/td corresponding to the diffusion time from con- 
tact A to B in our samples, and show the figures in table 
[J] Even though the sign reversal of the non-local con- 
ductance of sample III may be related to H/td, no clear 
features are seen at this energy for sample A and B. 

To conclude, we have presented non-local voltage 
measurements on superconductor-ferromagnet spin- valve 
structures, and observed a negative four-probe resistance 
that provides unambiguous evidence for crossed Andreev 
reflection as the dominating subgap transport process in 
some samples. This is an important prerequisite for the 
use of crossed Andreev reflection in efficient entanglers. 
The qualitatively different behavior reported by Russo 
et al. and also observed in one of our samples requires 
further systematic experimental investigation as well as 
theoretical work. 
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